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Abstract 

We derive the gravitational Hamiltonian starting from the Gauss-Bonnet 
action, keeping track of all surface terms. This is done using the language of 
orthonormal frames and forms to keep things as tidy as possible. The surface 
terms in the Hamiltonian give a remarkably simple expression for the total 
energy of a spacetime. This expression is consistent with energy expressions 
found in hep-th/0212292 . However, we can apply our results whatever the 
choice of background and whatever the symmetries of the spacetime. 
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1 Introduction 



The energy of gravitational systems has attracted much interest down the years (see, 
for example [ H| 2 ^Pl |i l |o 1 |o T| 7 || H lPllUpi 1|I12|I13|I14|I15| ) . In particular, for Einstein gravity 

S = kJ R-2A, (1) 

expressions were found for the energy of asymptotically flat pQ and asymptotically 
AdS [2 spacetimes. Some time later, Hawking and Horowitz [5] gave a general deriva- 
tion of the gravitational Hamiltonian, keeping careful track of surface terms. When 
evaluated on a solution, this Hamiltonian gave an expression for the total energy. 
This agreed with the earlier expressions found in [1,2 . However, this new expression 
could be applied regardless of the asymptotic behaviour of the solution. 

Recently, Deser and Tekin have found energy expressions for higher derivative 
gravities [T4*l ll5|. This inital work has focused on test spacetimes that are asymp- 
totically maximally symmetric, with background spacetimes (vacua), defined to be 
maximally symmetric everywhere. In this paper, we will perform the analogue of 
Hawking and Horowitz's calculation by deriving the Hamiltonian for Gauss-Bonnet 
(GB) gravity: a special, but important class of higher derivative gravities. We will 
obtain an expression for the energy by once again evaluating this Hamiltonian on a 
solution. As with j^j, our expression will be consistent with earlier results [T4*l ll5|. 
but can be easily applied even if the background is not maximally symmetric. This 
greater flexibility allows for a more natural choice of vacua in some cases, as we will 
illustrate with an example. 

Before going any further, let us briefly mention what GB gravity is, and why it is 
important. GB gravity is given by the addition of the Gauss-Bonnet invariant to the 
Einstein-Hilbert action 

S = k J R - 2A + a (R 2 - AR AB R AB + R A bcdR ABCD ) ■ (2) 

In D = 4 dimensions, the GB term is a topological invariant and does not enter the 
dynamics. This ceases to be the case in D > 4 dimensions. Furthermore, in Einstein 
gravity, the vacuum field equation are given by a linear combination of the Einstein 
tensor and the metric. In four dimensions, this is the most general combination of 
tensors that satisfies the following conditions [To] : 

• it is symmetric. 

• it depends only on the metric and its first two derivatives. 

• it has vanishing divergence. 

• it is linear in the second derivatives of the metric 1 . 
1 In D — 4, this condition is actually implied by the other three. 
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However, it has become unfashionable to restrict ourselves to just four dimensions. 
If we go to D = 5 or 6, it turns out that these conditions are satisfied by a linear 
combination of the metric, the Einstein tensor, and the Lovelock tensor fW\ \17\. The 
Lovelock tensor arises from the variation of the Gauss-Bonnet term in the above ac- 
tion (J2J). As D increases, we can use the same arguments to introduce higher order 
invariants (the Euler characters) [18J. In this sense, GB gravity is the natural gen- 
eralisation of Einstein gravity to 5 or 6 dimensions, although other higher derivative 
gravity theories have also been studied (see, for example [rmEU|01|.l22|. 0H] ) . 

It is interesting to note that the GB action (J2J) is renormalisable [211 , but perhaps 
the most compelling reason to study GB gravity is its appearance in String theory. 
Consider the slope (ie. a') expansion for the heterotic string. At lowest order, it 
is given by the Einstein-Hilbert action The first stringy corrections give rise 
to quadratic curvature terms (2211221 • For this effective theory to be ghost-free, the 
quadratic curvatures must be in the combination given in the GB action (j2J) [23 
I2H1- This link to String theory has generated a lot of research into higher derivative 

gravity and cosmology [21 E01EIlE21E3ISllE3EniEIlEHl 23110111111121113111113113 

iZlll3ll3inni- Although the early work j23E01EIlElEni achieved reduction to four 
dimensions via Kaluza-Klein compactifications, more recent studies |3~^l3*51l3~6U37tl38| 
IHHll4( )U4 !U42U4aU44U45U46U47U48U49U5 ( )j have focussed on the braneworld scenario 
Braneworld models are themselves inspired by String theory j2H| so it is natural 
to ask what effect any stringy corrections might have on their cosmology. From a 
holographic point of view, we might expect such higher curvature terms in the bulk 
to correspond to next to leading order corrections in the 1/N expansion of the CFT on 
the boundary /brane [13|. Calculating the GB Hamiltonian will allow us to investigate 
the GB version of "exact" braneworld holography . This will be discussed in 

a future article [55] . 

The rest of this paper is organised as follows: in section 2 we will give a deriva- 
tion of the GB Hamiltonian, with extra details given in the appendix. We will use 
orthonormal frames and differential forms to keep things as tidy as possible. Some 
readers may wish to ignore the details of this derivation and go directly to the energy 
expression at the end of the section. In section 3 we will show that this expression is 
consistent with existing literature (HZ1II11II3- We will illustrate the flexibility of our 
approach with a special example in section 13.21 Section 4 contains some concluding 
remarks. 

2 Derivation of the Gauss-Bonnet Hamiltonian 

2.1 The action 

The Gauss-Bonnet action (j2J) is most elegantly written in terms of differential forms. 
Suppose our D-dimensional spacetime, Ai, has metric, 

g = VAB E A ® E B (3) 

where {E A } is an orthonormal basis of 1-forms, and indices are raised/lowered using 
Vab = diag(— 1, +1, . . . , +1). We write {Xa} for the dual basis of vectors. 
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We will find it useful to define the following forms, 



e Al ...A m = , D ^ m y e Al ...A m A m+1 ...A D E A ™ +1 A ... A E A ° , (4) 

where eA 1 ...A D is the totally antisymmetric tensor with eo...(D-i) = 1- Notice that the 
scalar D-form, e, is the volume measure on Ai. 

Now suppose we have vanishing torsion, and that uj A b is the connection 1-form 
compatible with the metric g. The curvature 2- form is given by 

Q A B = du) A b + uj a c A uj c B = \r A bcdE c A E d . (5) 

The right hand equation above gives the Riemann tensor, R A bcd- The Ricci tensor 
is then defined by Rbd = R° bcd and the Ricci scalar by R = T] AB R AB - The GB 
action (J2J) can now be written 

S = k I -2Ae + tt AB A e AB + att AB A tt CD A e ABC D- (6) 
Jm 

where we have made use of the following identity [T%l l22|. 

E B A e Al ... Am = 5 Am e Al ... Am _ 1 - 5 B m _ l e Al ... Am ^ 2Am + ... + {-l) m ~ 1 5 Ai e A2 ... Am , (7) 

If Ai has a boundary, dAi, we need to define boundary conditions on dAi. We 
usually demand that the geometry of the boundary is fixed. However, as it stands, 
the action © is inconsistent with these boundary conditions. This is because its 
variation with respect to the metric does not vanish on shell. 

To cure this, we need to add a boundary integral. In Einstein gravity, this is 
the well known Gibbons- Hawking term jSHj. The generalisation to higher derivative 
gravities was worked out by Myers [THj. For GB gravity it is given by 



boundary = -« / AB A 6 AB + 2a0 AB A (V D - V ' E A ED ) A e ABC D, (8) 

JdM V 6 J 

where 9 AB is the second fundamental form [q^ IdTI]. It is defined as 

e AB = u AB - u AB . (9) 

where u AB is the connection for the product metric, go> that agrees with g on the 
boundary. The second fundamental form is closely related to the extrinsic curvature 
of DM in M [H]. 

The fully consistent action for Gauss-Bonnet gravity is therefore given by, 

S = 5bulk + ^boundary (10) 

where Shmik an d ^boundary are given by equations © and (JBJ) respectively. It is inter- 
esting to note that the boundary integrand contains a linear combination of the first 
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two Chern-Simons forms. This is related to the fact that the bulk integrand contains 
the same linear combination of the first two Euler characters [T5] . 

Finally, we note that although the action (jlUj) is well defined for spatially compact 
manifolds, it is divergent when Ai is spatially noncompact [3] . To get round this we 
need to choose a reference background, M. with metric g. This background should 
be a static solution to the field equations = 5J, but does not have to be maximally 
symmetric. The boundary conditions are unchanged which means that dAi should 
have the same geometry as dAi. We can then define the background action, S, in 
the same way as before. The physical action is the difference, 

I = AS = S - S. (11) 

In the Hamiltonian picture, this background can be thought of as defining a back- 
ground energy or zero energy solution. For example, for an asymptotically AdS 
spacetime, we would probably choose the background to be pure AdS, but we do not 
have to. Any asymptotically AdS black hole spacetime would be equally valid. 

2.2 Splitting space and time 

Although the physical Hamiltonian will be constructed from the action (fTTj) . it is 
clear that it will just be the difference of the Hamiltonian constructed from S and 
that constructed from S [Hj. For the time being we will concentrate on the former. 

2.2.1 Foliations of M 

To proceed, we need to deconstruct the spacetime Ai by separating space from time 
in the following way. First, we choose a timelike vector field, d/dt. Now introduce 
a family of spacelike hypersurfaces {£t} labelled by the parameter t. This family is 
a foliation of the full spacetime. We assume that the hypersurfaces have no inner 
boundaries and do not intersect each other. They meet the timelike part of the 
boundary (call this B) orthogonally, and in the far past/future, they coincide with 
the spacelike part of the boundary (call this £oo)- Therefore the total boundary, 
dM = B U Soo. 

We can write the metric for Ai in ADM form pQ, 

g = -N 2 dt 2 + 7a6 (t, x a )(dx a + N a dt)(dx b + N b dt), (12) 

where N is the lapse function, N a the shift vector, and / j a b{t,x a ) the induced metric 
on E t . It is natural to choose the following ortho normal basis of 1-forms, 

E L = Ndt, E a = E a b (dx b + N b dt) (13) 

where 5 a bE a c E b d = •y cd . We would like to emphasize some notation at this point. 
Lowercase latin indices label components in E t , whereas uppercase latin indices label 
components in M.. For example, {E a } is an orthonormal basis for T, t whereas {E A } = 
E 1 - U {E a } is an orthonormal basis for Ai. 
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The dual basis of vectors is given by [OT] 



1 (1 - N"— 

N \dt dx° 



X„ = E, 



d 



a dx a ' 



(14) 



X± is the vector normal to S t , and is not necessarily equal, or even parallel to d/dt. 
We should also note that X± is tangent to B, as T, t and B are orthogonal. 

Normally, the next step is to use the Gauss-Codazzi equations jH2] to rewrite 
the bulk part of the action (see for example [3]). This has infact been done for GB 
gravity (03] although the contribution from surface terms was ignored. In this paper 
we are using the language of orthonormal frames and differential forms. We therefore 
need to know how to translate the Gauss-Codazzi equations into this language. This 
is explained in [64] . so we will merely state their results. 

The Gauss-Codazzi equations describe the decomposition of the bulk Riemann 
tensor into spatial tensors defined on H t - In the language of forms, it is the curvature 
2-form that we wish to decompose. We start by decomposing the connection, uj ab . 



to 



.-La 
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lib 



-H a + a a E 
u ab + l ab E~ 



(15) 
(16) 



Here, we have two 0-forms: a vector, a a , and an antisymmetric tensor, l ab . The 1- 
form H a = H a i ) E b , where H a b is the extrinsic curvature of H t in Ai. From now on 
anything labelled with a tilde is intrinsic to as opposed to M.. Therefore, Co 
the connection for the induced metric 7 = S a bE a ® E b . 

We now write the Gauss-Codazzi equations in the following way 



~ab ig 



Q La = -VH a + E L A 
n ab = Q<*> + H a AH b + E ± A 

j~ab 



-$LH a - ^V(iVa a ) - l ab H b 



1 ~ 



ii^ - — V(Nl ab ) + H a a b - H b a a 
N 



(17) 
(18) 



where Vt ab = du ab + u a c A u cb is the curvature 2-form for E t , and the operator V is 
the covariant exterior derivative [59(150] (on E t ). Note that V reduces to the exterior 
derivative, d, when acting on scalars. The definition of the operator $^ is given in j64j . 
It is closely related to the Lie derivative with respect to the vector X±, although it 
lives entirely on S t and acts on tensor components as if they were scalars. In many 
ways it behaves like a partial derivative 2 . 

In principle we could also decompose the torsion 2-form (see JOl]). However, we 
have set torsion to zero, which means that every component of the torsion decompo- 
sition must vanish. This gives the following conditions: 



a a E a 

VE a 

$ ± E a 



H ba , 

dN 

~N ' 




-H a - r h E b . 



(19) 

(20) 

(21) 
(22) 



It is often useful to think of $ a as the frame-form version of 

at "t 
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2.2.2 Foliations of B 



Since we intend to keep careful track of surface terms, we will need a foliation of B, 
as well as M.. On B, the foliation is given by the family of surfaces {St}. For a given 
value of t, St is the intersection of B and E t . 

We need to understand how bulk quantities project on to S t . Near St, the metric 
can be written in Gaussian normal coordinates 

g = dz 2 - E 1 ® E L + S ij E i <g> BP (23) 

From now on we will write E z for the extension of dz into the bulk, and X z for the 
extension of the inward pointing normal d/dz. Notice that we are using indices i, j 
etc to label components in St- 

We can use the techniques developed in jM] to project bulk quantities on to B 
and then on to St- We find that the decomposition of the connection is given by 

(24) 

iE z (25) 
(26) 

C j E z (27) 

where \^ is some antisymmetric tensor, and 

c = a z (28) 
h l = l zi = —H zi (29) 

Note that anything wearing a hat is intrinsic to S t - Therefore, cD*- 7 is the connection 
for the induced metric (on S t ), A = S^E 1 ® EK We have also defined, 

ir = IPjBP, (30) 
K l = K^EK (31) 

Here, we should interpret H % j and K % j as the extrinsic curvatures of St in B and S t 
respectively 



u ±z = 


hE l 


+ 


cE 1 - 




u ±l = 




+ 


a i E ± 


+ 


u zl = 


k l 


+ 


VE 1 




J j = 




+ 


l ij E^ 


+ 
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The curvature decomposition is given by 

n Lz = d^E 1 ) +H* AKi + E z A{---} 



+E 1 - A 



1 



$1 {hE l ) - -d(Nc) - a l Ki - b'H, 



n- 



n z 



-VH l + bjE j A JC + E z A {• • ■ } 



+E L A 



hH* - — V(AV) - l ij Hj + dC 
N 



UbjE 1 



VK i -b j E j AH i + E z A{---} 



+E A 



_iT - — V(AV) + K/ j - fc^V - cff 1 ' 



+ if* A AF - AT* A K j + E z A {• • ■ } 



+E A 



S-lw" - jjV(Nl ij ) + 2H*a? ] + 2K [i V ] 



(32) 



(33) 



(34) 



(35) 



where we will not care what is contained in {■ ■ • } as the integration of E z over B is 
zero. In analogy with the previous section, Q 1 ^ and V are the curvature form and 
covariant exterior derivative on St, respectively. Again, V reduces to the exterior 
derivative, d, when acting on scalars. The operator $^ is the analogue of $j_ on B 3 . 

As long as we are near B, equations (|3*2*|) to (|33|) are the frame-form version of the 
Gauss-Codazzi equations for a hypersurface of codimension two 



2.3 The Hamiltonian 

We are now ready to start calculating the Hamiltonian. However, we will find it 
convenient to continue working with the action until virtually the bitter end. When 
our action finally has the desired form we will switch to the Hamiltonian picture, and 
give an expression for the gravitational energy of a solution. 

We will start with the bulk part of the action (jBJ). Our aim is to write it so that it 
contains no derivatives of the lapse function or the shift vector. This is because these 
are ignorable coordinates, and should behave like Lagrange multipliers. They will be 
paired with the Hamiltonian and momentum constraints respectively, as is the case 
in Einstein gravity [3]. We also want to eliminate terms like $j_H a , which contain 
second time derivatives of the canonical variable E a . We will need to use integration 
by parts to achieve these aims. This means that the bulk action (JHJ) will contribute 
surface terms. In summary, we expect to write © as 



5bulk — 'S'bulk ^leftover; 

where SWtover are the leftover surface terms, and 



°bulk 



K 



dt 



7r a A E a — NtC — N a H a . 



(36) 



(37) 



3 If A is an arbitrary p-form in S t , then near St , we can write A = A + A z A E z , where A and A z 
arep and (p— l)-forms respectively, living on St- It can then be shown that §±A = %±A + E Z A{- • • }. 
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Here n a is the momentum conjugate to E a , and 7i and H. a are the Hamiltonian and 
momentum constraints respectively. We have also introduced the intuitive notation 

E a = $a_E a . 

at 

At this stage, the surface part of the action is given by equation (jHJ). We can split 
this into two parts, 

"^boundary = Sb + "Soo, (38) 

where Sb contains the integration over B and the integration over E^. In sec- 
tion ^. 51 we will group these terms with SWtover- This will give us a modified boundary 
term 



q* — Q A_ Q 

boundary "-"boundary T "-"leftover- 

which will be closely related to the gravitational energy of a solution. 



(39) 



2.4 The bulk 



As promised, we start with the bulk part of the action (jUJ), in the hope of deriving 

5 buik and 
given by 



S^ ulk and S'ieftover- Making use of the formula (J7J), we find that the bulk action © is 



where 



S 



bulk 



^kinetic + 5*1 + S2 + S3 



(40) 



with 



^kinetic = K I E f\ \ —2$±H A (b 
J M 



-4a 



VH b A $ ± u cd + $ ± H b A F 



b a zpcd 



AC 



bed 



Si = k J E^ A I -2AC + F ab A U + aF a0 A F ca A ( abcd 

2 ~ 



S 2 



M 



k E A < -— V(Na b ) A C 



N 



-4a 



1 ji 



N 



V(Na b ) A F cd + 2VH b A H [c a d] 



(41) 



(42) 



k E A < 4a 
'm 



VH b A ^V(M cd ) - l ba H a A F cd 



AO 



A (bed } (43) 
(44) 



bed 



na6 oafe 



ft a6 + H a AH b 



(45) 



Note that we have introduced the analogue of eA 1 ...A m on S t , 

1 



Cai- 



■ On 



(D - l-m)\ 



£a 1 ...a m a m+1 ...a D _ 1 E ai A ... A E" 10 ' 1 . (46) 
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We will now rewrite expressions (jlTJ) to (|44|). bearing in mind the goals we mentioned 
at the beginning of this section. The kinetic term is given by 



5 k i nc tic = « J dt J n a AE a - N a H a 

+k J dt J AaNH b A $ ± u cd A (bed + (-l) D N a n b E b a - S^. (47) 

where 

7i a = -2( ab A H b - 4a( abcd A H b A A cd (48) 

A cd = Q cd + -H c AH d (49) 
3 

and the Hamiltonian constraint, 

W a = (-I^-^.Vtt*. (50) 

In section we will show that ir a is indeed the momentum conjugate to E a . The 
derivation of equation (J47|) is given in the appendix. 

The remaining terms in the bulk action can be written as follows: 

Si = — k f dt f NH (51) 



Hi 



S 2 = k J dt J N [2a b Cb + 4:aa b F cd A Q bcd ] 



(52) 



S 3 = — k I dt I AaNVH b l cd A ( bcd (53) 



bicd 
<±u:i v v J-i I 

St 

where the Hamiltonian constraint, TC, is given by 

H = 2A( - F ab A ( ab - aF ab A F cd A ( abcd . (54) 
Note that in deriving the expression ()53|) for 53, we made use of the following: 

V 2 H b = Q b c A H c . (55) 



H n A 



A ( bcd = 0. (56) 



Collecting together equations (|T7j) . (joTj) . ()52jl. and we see that Sb u ik takes the 
desired form 

'S'bulk = 'S'bulk ^leftover (57) 

with S£ ulk given by equation (}3*T|) and SWtover given by 

leftover = -S»+ K J dt J (-1)° N a 7T b E b a 



+ N^2a b ( b + 4a H b A $±u cd + a b F cd — VH b l cd A ( bcd ^ . (58) 
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2.5 The boundary 

As expected, rewriting the bulk part of the action in the desired form (J37|) has altered 
the boundary part of the action. In particular, we have a leftover surface integral (|58|) 
that must be added to the original boundary part of the action (JHHj) . This gives the 
modified boundary action, 

'-'boundary = ^leftover + Sqq + Sb- (59) 

In order to combine each term in the above equation, we need to write them in a 
common form. This will involve integrations over S t , of well defined quantities on St- 
Let us begin with Stover + Sao- From we know that near B, E z a = 5 z a and 
N z = 0. This means that the momentum term in (|58|) gives 

k dt (-l) D N a n b E b a =k I dt I (-l^JVVj-ETV (60) 
J J s t J J St 

We now use the fact that 

%jo zi = % ± k' + E Z {■■■}. (61) 
along with (J28|) and (|29|) to rewrite the remaining terms. We find that 



leftover + S^ = K \ dt / (-l)^iVV^ + 2iV C 



+ A01N 



St 

1 A4>ij (62) 



-b k E k A + cF 11 + F z l 13 - 2H l A $ ± K 3 - 2a l F z > - 2F t V 



where F b = —VH b . We have also introduced the S t analogue of C<n-a n an d &Ax-A n - 

Meanwhile, terms like F z and i™ can be deduced by comparing equations (|32j) to 
fl35j) with (JT7D and (HH1). 

F z = d (b k E k ) + H k A K k (64) 
F i = -VH* + b k E k AK* (65) 
F z i = VK j -b k E k AH j (66) 
pij = ffi+^AF-^AF (67) 

Note that we have dropped all terms like E z {- • ■ } as we are now integrating over B. 
Now consider Sb- The only non-zero components of 6 AB on B are: 

e ±z = to ±z , e zi = u zl (68) 



n 



where to and uo zl are given by equations ()24|) and ()26|) respectively. We use this 
fact, along with equations and (fHo"|). to write Sb in the following way: 



S B = k I dt I 2NK 1 A (pi- 2iVc0 + AaNK 1 A 

2 

% 



3 



A 



+4aiV ^ -26 



F 3 - -&,£ fc A 



— c 



F 1 - 7 + -7T A K 3 
3 



+ 27P A 
+ & fc £ fe A 



G J + - ( -ctf 3 ' + &^,£ fc 



+ -b [i k j] 

3 



- — V(W) - P k H k + cK j - Vb k E k 



(69) 

o j l ° J ) 

where 

(70) 
(71) 

The expressions (j62j) and (J6"9"|) now have a common form, so we can combine them to 
get 

(72) 



G ij = $ ± Q ij - — V(Nl lj ) + 2H [i a j] + 2K [i b> ] 
N 



'-'boundary ^boundary,! ^boundary, 2 ^boundary, 3 



where 

'S'boundary,! 



k I dt I (-l) D N l 7TiE 3 i + N i 2K l A 



3^ » 



St 



+AaK l A 



1 



Qjk + H j A R k _ _ Rj A K k 
3 



A (j)ij k \ (73) 



S ] 



boundary, 2 



J dt j AaN jd (b k E k ) l ij - -jyhE k A S/(NP 



Si 



boundary, 3 



-2a l VK 3 - —K l A V(AV) \ A r „ 
k I dt I AaN \ -2H* A $ , K j - 2K l A $ , H j 



+H k A - 2K l A F fc F fe S> A <b u 



Since we have zero torsion, we use the condition 



V0i r .. in = 



(74) 



(75) 



(76) 



to show that the integrand in S DO undary,2 is a total derivative. Since St is a boundary 
of Et, it has no boundary of its own. This means that, 

^boundary, 2 = (77) 

By integrating by parts on the time derivatives in (fT5J) . we can also show that 

S b oundary,3 = « / dt j 4«iV |-2K l A W A H k A . (78) 
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Here we have used the zero torsion result 

$ ± E { = -H* - VjEP. (79) 

and the identity 

H k A K k l lj A <j>ij - 2K l A H k l jk A - 2K l A H j l\ A E l A <j) ijk 

= K i AE i A H j l kl A (f) jkl = 0. (80) 

The last equality in ()8U|) follows from the symmetry of K l A 

We now collect together equations (J75|) . (J77)l and (|7Hj) to deduce that, 

boundary = « / ^ ^ (-1)^7^ + N [2K* A fa 

' ~' A&jfcl (81) 



+4aK i A 



fij* - #J A if fc — -K j A K k 
3 



2.6 The conjugate momentum 

Now that we have the action in its correct form, it remains to calculate the momentum 
conjugate to E a . This is given by, 

p - = (82) 

where £buik is the bulk integrand. We could take £buik from S^ ulk . However, it is 
convenient to temporarily undo the integration by parts that gives equation (|148|) 
(see appendix). In other words, we leave the derivatives of N a in the bulk action. 
This is perfectly OK, because it does not affect the bulk dynamics, and therefore the 
value of the conjugate momentum. The bulk integrand is temporarily given by, 

£bulk = TTa A N$ ± E a - NH (83) 

where 7r a and H are given by equations (|4*Hj) and ([54)) respectively. Using the zero 
torsion decomposition we can say 



^buik = -Nn a A (H a + l a b E b ) -NH= —Nit a A H a — NH (84) 



where the right hand equation follows from ([149)1 . and the antisymmetry of l ab . 
Referring to equations ()22|) and ([146)1 . we use the chain rule to show that, 



1 ^bulk /q r\ 

Pa = ~K^fT~ = 11(1 

N dH a 



This non-trivial result is due to the following cancelation, 



A H + = 0. 86 
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We conclude that ir a is indeed the conjugate momentum. It should be thought of as a 
function of E a . To derive the Hamiltonian, we should invert this function. However, 
7r a is cubic in E a , so the inverse is multivalued. This is a well known property of higher 
derivative gravities. In the Hamiltonian picture, this could mean that we could jump 
from one solution to another. These "zigzagging" histories still provide an extremum 
of the action. In this paper, we will assume that at any given time, we have a unique 
solution. This is just the same as saying that we are not in the process of jumping 
from one solution to another. For more discussion on multivalued Hamiltonians in 
this context, refer to 



2.7 The physical Hamiltonian 

We have shown that we can write our action as, 



$ — ^bulk + ^boundary — K I dt 



where 

C* hnlk = n a A E a - NH - N a H a 
N 1 2K { A <f)i + AaK' b A 

+ (-l) D jV i 7T j £' i . 

The corresponding Hamiltonian is defined as, 
H 



^bulk 



-^boundary 



C 



boundary 



1 



Qjk _ R j A R k _ _ R j A R k 
3 



A0, 



(87) 



(88) 



19) 



K TT a A E a £ bu i k K / ^boundary 

= k [ NH + N a H a - k [ £ boundary (90) 

JT. t J St 

To arrive at the physical Hamiltonian, we need to subtract off the background Hamil- 
tonian, 



H 



k / £ 



boundary 



(91) 



Here we have used the fact that the background is a stationary solution to the field 
equations 

H = Ha = Tta = 0. (92) 

The physical Hamiltonian is therefore given by, 



« / NH + N a H a -K / AC 



phys — n I Jv/ltiv /L a n / ^-^boundary 
St J St 



(93) 



where 



A^Lundary = ("1)^7^ + N 2AK* A 0, 



+4a 



AK { A (h jk - H j A H k ^j — -A (R* A K j A K k 



A<f>ijk\ (94) 
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For any quantity Q in the test spacetime with corresponding quantity Q in the back- 
ground, AQ = Q — Q. Notice that we have Af2 jfc = AiP = 0. This is because the 
geometry of the boundary is the same in the test spacetime, as in the background. 
If our test spacetime is a solution to the field equations, it satisfies the constraints 



n = n a = o 

Its energy is then given by the value of the physical Hamiltonian, 



(95) 



E 



-K 



1)'\Y '-,/•;', + N { 2AK l A 



s t 



+4a 



A/v ' A ( W k - H j A H k ) — - A [K 1 A K j A K k 



1 



A 4>ijk \ (96) 



Given the technical complexity of Gauss-Bonnet gravity, we believe that this expres- 
sion is remarkably simple. Note that for a = 0, we recover the correct result for 
Einstein gravity, as of course we should. 



2.8 Using a coordinate basis 

Although the final result (|9T)j) is neat and tidy, we might prefer to work in a coordinate 
basis, and express the Hamiltonian in terms of the familiar tensors of General Rela- 
tivity. In this case, our canonical variable is the induced metric 7 a fc. The conjugate 
momentum, TT ab , is given by j^l], 

With this in mind, we can verify that the Hamiltonian and momentum constraints 
given by equations (J53)l and ([5U]) respectively, agree with the corresponding expres- 
sions in [63J. If z A is the normal to the timelike boundary, B, and n A the normal to 
S t , the Hamiltonian (}9*B]) can be written 

= K / d D ~ x x [NH + N a H a ] - k [ d D ~ 2 x AC* houndaiJ (98) 

JT. t JSt 

where now 
A^bounda ry = 2N%z j + NV\{2AK 



AKi (R? k mn - 2H^H k ) - ^A (Jf/Jf^) 



(99) 



Here, B? k m n is the Riemann tensor on S t , constructed out of the induced metric A^. 
Kj and Hj are the extrinsic curvatures of St in £$ and B respectively, and K is the 
trace of Kj . More precisely, 

Kn = -X^X^Zm, H id = -X%X^Vm m , K = X\K). (100) 
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3 Comparison with previous energy expressions 
3.1 Deser and Tekin 

Now that we have derived an expression for the energy in GB gravity, we should com- 
pare it with previous results in the literature. In particular, Deser and Tekin [T4*l lTo] 
used a "conserved charge" technique to derive the energy of asymptotically maximally 
symmetric spacetimes above backgrounds of constant curvature. This method can be 
applied to generic higher derivative gravities, but we can check it is consistent with 
our result in the Gauss-Bonnet case. 

Suppose our test spacetime Ai, is asymptotically maximally symmetric. We 
choose our background, Ai to be the maximally symmetric solution with curvature 
form, 



The effective cosmological constant is given by, 

1 



A. 



eff 



1 ± Vl + 8Aa* 



„ (£>-3)(P-4) 

a = a-— — r. 102 

(D — 1)(D — 2) v ; 



which is real and multivalued for 1 + 8Aa* > 0. 

Our aim is to calculate the energy associated with the timelike Killing vector 
d/dt. We can choose our foliation {E t } so that the shift vector vanishes. Since d/dt 
is Killing, it is clear that H l = 0. The expression for the energy reduces to, 



h = -ii J Ni2AK' Ao,+ in 



AK { A h jk A I K { A K j A K k 

3 



A <pij k 



(103) 

Since S t lies in the asymptotic region, we can say, 

A (k* A K j A K k ^j A (f) ijk w 3AA^ A (R j A K k ) A (f> ijk , (104) 

where Ri is Ri evaluated on the background. The energy is now given by, 

E = -kJ N ^2AK { A <j>i + AaAK* A (tt jk - K j A K k ^j A ijfc } (105) 

Using the fact that H 3 = 0, it is clear from equations ()35|) and (|101j) that, 

2 A 

W k - & A K k = {D _ 1) f D _ 2 f j A E k . (106) 

This implies that the energy, 

E = —k (1 + 4a*A c fj) / 2NAK i A fa = ±k<A + 8Aa* / 2NAK i A(p i (107) 

J s t J s t 

where we have used the S t analogue of the useful formula (|7ji. and the cosmological 
constant relation (|102|) . 
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In order to make contact with [T^ ITo]. we switch to a coordinate basis, 

E = ±«Vl + 8Aa* / d D ~ 2 xV\2NAK (108) 
J s t 

We now follow the procedure described in |3] for Einstein gravity. Let us start with 
the test spacetime. Near St, we can express the metric on S t in Gaussian normal 
coordinates, 

ds| t = ^ ab dx a dx b = dz 2 + q^z, x k )dx i dx j . (109) 

where (^(O, x k ) = Xy is the metric on St- Similarly, for the background, we can write 
the metric on S t (near S t ), as, 

ds| = 'j a bdx a dx b = dz 2 + q~ij(z,x ■ )dx l dx J . (HO) 

To ensure that the normals to St agree on the test spacetime and the background, we 
choose the diffeomorphism z = z and x 1 = x l . In these coordinates, 

K=-~qVq ij>z , K = ~<fiq ij;Z . (Ill) 

Since both metrics agree on the boundary, we note that Aq^ = there. Therefore 
on S t , 

A^ = -^(A ? ,)^ = ^(A 9 ) i2 . (112) 
where Aq = X^Aq^. This gives a final energy expression, 

E = tkVI + 8Aa* / d D ~ 2 xVXN ( Aq) z . (113) 
J s t 

We now use Deser and Tekin's method [r^l lTo"] to calculate 4 the energy associated 
with the timelike Killing vector, t A . 

E DT = T kVI + 8Aa* J dS A n B {t c V B h AC - t c V A h BC + t B V A h - t A V B h 

+h BC V A t c - h AC V B t c + t A V c h BC - t B V c h AC + hV B t A ] (114) 

where h AB = g AB — g AB and h = g~ABh AB . Here qab and qab are the metrics on 
M. and M. respectively. We will choose to work in a synchronous gauge for which 
n A hAB = jSJ. As the metrics agree on the boundary we can also set h AB = on S t . 
If St has inward pointing normal d/dz, the measure is given by dS A = -d D - 2 xVX5 z A . 
For vanishing shift vector, we have t A = Nn A , and the energy ()114|) simplifies to, 

E DT = ±kV1 + 8Aa* [ d D ~ 2 xVXN (d b h zb - d z h) . (115) 
J s t 

4 In |14II15| . the authors do not explicitly write down an energy expression for the GB action with 
a bare cosmological constant A. However, they give enough information to easily derive equation 

033}. 
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In the Gaussian Normal coordinates we have recently described, h zb = and h = Aq. 
Deser and Tekin's energy now reads, 

E DT = =fm/1 + 8Aa* [ d^xVXNd 2 (Aq) . (116) 

J s t 

This expression clearly agrees with the equation (|113|) . We can conclude that although 
our derivation was very different to that in OHH], our results are consistent. 

3.1.1 Application to GB black holes 

One of the nice features of GB gravity (J2J) is that it contains static, spherically sym- 
metric solutions jnaEainHiinniizniiziiiiniinni of the form, 

ds 2 = -V(r)dt 2 + y^r + r 2 dn 2 D _ 2 (117) 

where dfl 2 D _ 2 is the metric on a unit (D — 2)-sphere. We will assume that l+8Aa* > 0, 
so that there are two possible branches for the potential, 



Here /x > is a constant of integration that gives mass to the spacetime. The upper 
branch has a naked singularity at r = 0, whereas the lower branch is a real black 
hole with a unique event horizon surrounding the singularity jSZl- We wish to cal- 
culate the mass, M, of these spacetimes above the appropriate maximally symmetric 
background, 

— dv 2 

ds 2 = -V(r)dt 2 + T — + r 2 dn 2 D ^ 2 (119) 
V(r) 

where 

v ^ = 1+ 2(d- m> - 2 )a . i 1 ± ^™^) ■ (120) 

The foliation in each case is given by surfaces of constant t, so we soon see that, 

Qij = r2 ( z )Xij, Qij = r 2 (z)xij (121) 
where Xij is the metric on the unit (D — 2)-sphere, and 

If r (0) = f (0) = R, on S t , it follows that, 

(Aq) tZ = - 2(jP ' " 2) A^/V(R) » -£zjW(J2) on S t . (123) 
R RVV 
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For large R, 



Finally, we note that N 



AV(R) 



1 



~Vl + 8Aa*R D - 3 
V and apply equation (|113j) to give, 

M = nVL D _ 2 {D -2)h 



(124) 



(125) 



where is the volume of the unit (D — 2)-sphere. This is the standard result. 

It is always valid for A c g < 0. For A e g > our analysis is valid only if the de Sitter 
horizon is much larger than the black hole horizon [To] . 

3.2 A special case 

In the last section, we assumed that 1 + 8Aa* > 0. Now consider what happens 
when 1 + 8Aa* = 0. We cannot make use of the expression ()113|) because it involves 
multiplying an infinite integral, by zero! We will not worry about how one would 
modify the approach of [T4llTo] to accomodate this. Instead, we will sell the approach 
developed in this paper. Let us focus on the 5-dimensional black hole with, 



V{r) 




(126) 



To calculate the mass, we need to go back to equation (|9H|) . We also need to choose 
a background. In this example, the maximally symmetric solution with \i = is 
not necessarily the most natural choice. We might prefer /i to be chosen so that the 
horizon has zero area jHZj. Whatever our choice, we illustrate the flexibility of this 
work by allowing for non-maximally symmetric backgrounds. We will keeps things 
general and say that the background potential is given by 



V(r) = 1 + — - 
4a 




(127) 



As before, our foliation is made up of surfaces of constant t, with boundary St given 
by r = R. To apply the energy expression (|96p. we need the following ingredients, 



N 



V, N l = 0, K l 



R 



E' 



E i = 0, VL jk = —E j A E k . (12* 
R 1 



Now use the useful formula J7J) in ()96|). to derive the energy, 

E = -Kn 3 R*Vv^A(Vv) + ^- A(Vv)-±A(yVv) } (129) 

To keep things tidy, we write 

V = y 2 - m, V = y 2 -m (130) 
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where 




m = 




m = 




(131) 



Now for large y, 




) 



(134) 



(132) 



(133) 



If we plug this back into ()129|) we find, 



E = 3kQ 3 [2aA(m 2 )] + O 




) 



(135) 



Now we send R, or equivalently y, to infinity, to derive the black hole mass, 



If we had chosen the background to be the black hole of zero size, we would have 
/2 = 2a. Our black hole mass would be given by M = 3kQ 3 (h — 2a), which agrees 
with the "minisuperspace" method employed in |57j . 

4 Discussion 

In this paper, we have derived a neat and easy to use expression for the gravitational 
energy of a solution in Gauss-Bonnet gravity. This was done using a Hamiltonian 
approach, much like the one used by Hawking and Horowitz jH] for Einstein grav- 
ity. Given the technical complexity of the derivation, our final expression is 
remarkably simple. 

There have been other ways of calculating the energy of certain Gauss-Bonnet 
solutions [HZIE1E3- We have shown that our Hamitonian approach yields results 
that are consistent with these. Each approach has its advantages and disadvantages, 
as we will now discuss. 

Consider the "conserved charge" method given in [TH ll5j. The authors identify 
a conserved current associated with a timelike Killing vector. The gravitational en- 
ergy corresponds to the "charge" of this current. This method can be applied to 
generic higher derivative gravities, of which Gauss-Bonnet gravity is just a special 
case. However, the background spacetimes, or vacua, are always assumed to be maxi- 
mally symmetric everywhere. That is not to say that this method cannot be extended 
to a more general choice of background. This should clearly be a topic for future re- 
search. It would also be interesting to know how to apply this method to the special 
case discussed in section I3~21 



M = 3kQ 3 A//. 



(136) 
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Similarly, we should also ask if we can extend our Hamiltonian approach to more 
general higher derivative gravities. It should be fairly easy to consider the Lovelock 
action [To] . 

m , 

S = V a n S n , S n = / A ... A Q A " B - A e AlBl ... AnBn . (137) 

n=o J M 

Here we have a linear combination of the first Euler characters, with the surface 

terms given by the same combination of Chern-Simons forms JH] . Life would be more 
difficult if we wanted to consider an arbitrary combination of Riemann tensors, as 
the surface terms are generally unknown 5 . 

The "minisuperspace" method used in [2Zj is closest in spirit to the Hamiltonian 
approach. The idea is to consider a static, spherically symmetric ansatz for the metric, 
and insert it back into the action. The action becomes one-dimensional, making it 
easier to fix the boundary term. When we turn to the Hamiltonian, and evaluate it 
on one of the black hole solutions given in sections 13.1.11 and 13.21 we derive the black 
hole mass. This method is very simple and easy to use, but somewhat limited. It can 
only be applied when the one-dimensional "minisuperspace" model is valid. This is 
OK for the black hole spacetimes discussed in [57], but a more general approach is 
clearly desirable. 

The Hamiltonian approach developed in this paper is the appropriate generalisa- 
tion. It can be applied whatever the symmetries of the solution, and without having 
to reduce the number of dimensions. In particular, we will use it to investigate the 
generalisation of braneworld holography [HE] for Gauss-Bonnet gravity [56J. 

To sum up, we have derived an expression ([Ho*]) for the energy of a solution to 
Gauss-Bonnet gravity. This can be applied whatever our choice of background, and 
whatever the symmetries of our solution. This should, hopefully, give us a platform 
to investigate Gauss-Bonnet gravity more thoroughly. 
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A Rewriting Skinetic 

We begin with the expression for Emetic given in equation (|41j). This contains terms 
like $±H a , which need to be eliminated by integration by parts. Using the fact that 

$±U... a m = $±£ am+1 A ( ai ... am+1 , (138) 

5 Attempts have been made to derive surface terms for general higher derivative gravities using 
auxiliary fields |?21EI1E]- However, in the GB case, these results do not agree with ^S], so one 
should proceed with caution. A deeper understanding of boundary terms versus boundary conditions 
can be found in |75j . 
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we find 

^kinetic 



k / E ± A { n a A $ ± £ a + 4a 



# e A $ ± tt cd - VH b A $ ± u 



~,cd 



AO 



/, / 2H b A Cb + 4aH b A A cd A Ctcd 



(139) 



where 7r a is given by equation (|4*KJ) and A cd by equation (jj§|) . 

If we note that 9 ±b = —H b is the only non-zero component of 9 AB on Eqo, we can 
use equations (jHJ) and (fTSj) to show that, 



5, 



x k I 2H b A ( b + AaH b A A cd A ( bcd . 



(140) 



This will cancel off the second line in equation ()139j) . 
We will soon need the following identities, 



V( av .. am EE EE VQ Cd . 



(141) 



The left hand side follows automatically from the zero-torsion condition (}2~T]) . whereas 
the right hand side is just the frame-form version of the Bianchi identity. We now 
make use of ()141|) . and the relation, 



to show that, 



E 1 A 



' M 



E 1 - A $ L Q cd = dtAV (N$ ± u cd ) 



A (bed 



(142) 



H b A $ ± & d - VH b A $ ± cD 

= J dt J H b A V (N$ L u cd ) - VH b A N% 
= - J dt J d [NH b A % 1 _Q cd A Cbcd] 
= J dt j s NH b A $j_uj cd A Cbcd 



~ cd 



(143) 

Note that we have applied Stokes' Theorem in the following way. If A is a (D— 2)-form 
on H t , then 



eft / dA= dt AdA = dt AdA 



- / d(dtAA) 



dtAA = - dt A (144) 



where we have used the fact that dA = dA + • • } jM]. We now insert (|14U|) and 
(|143j) into equation (|139|) to give, 



S ] 



kinetic 



k [ dt [ n a A N$±E a + k f dt f AaNH b A$ ± Cu cd A ( bcd - S^. (145) 
J Jz t J J s t 
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Notice that $±uj cd has been removed from the bulk part of Skinetic- This is because 
we have zero torsion. u cd is not an independant dynamical variable, so any time 
derivatives of it should indeed disappear from the bulk. 

We are not yet finished with Skinetic- Equation (jl45j) still contains derivatives of 
the shift vector. From equation (JHJ), we deduce that 

$ ± E a = — (E a -$^E a ), N = N a — = E a b N b X a . (146) 
N \ J ox a 

Since N lives entirely on Sj, $^ is just the intrinsic Lie derivative on S t [S3]. There- 
fore, 

%$E a = i${dE a ) + d(i^E a ) = -N% d i {w\)E c + V{E a b N b ). (147) 

where iyA is the interior 'product of the vector Y and the p-form A 76 . After some 
integration by parts we see that, 

dt [ ir a A$^E a = [ dt [ N a H a - f dt f {-l) D N a n b E b a , (148) 
is t J is t J Js t 

where the momentum constraint, 7i a , is given by equation (jSUJ). In deriving this 
constraint, we have used the fact that, 

7T b A E c = 7T C A E b . (149) 

This is not obvious but can be shown using the symmetries of H a b and the Riemann 
tensor 6 (on S t ). 

By inserting equation (j!48)l into ()145j) . we arrive at the final expression (jlTjl for 

Skinetic • 
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